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Abstract 

Current  guidance  and  control  systems  for 
interceptor  missiles  use  algorithms  that  require 
knowledge  of  the  target.  These  systems  are  not  robust 
to  target  maneuvers  and  disturbances.  In  this  paper,  we 
investigate  the  use  of  a  novel  controller  named  a 
Sliding  Mode  Controller  (SMC).  We  consider  two 
types  of  interceptors  or  kinetic  energy  kill  vehicles 
(KEKVs):  a  ballistic  missile  and  a  cruise  missile. 
Sliding  mode  controllers  are  robust,  simple  guidance 
and  control  algorithms  which  can  be  applied  to 
interceptors.  After  the  interceptor  flight  path  angle 
captures  the  line-of-sight  angle,  the  SMC  will  follow 
the  line-of-site  angle  regardless  of  external 
disturbances,  plant  uncertainties  or  measurement  noise 
until  the  target  missile  is  intercepted.  The  developed 
guidance  and  control  algorithms  require  estimating  or 
measuring  only  the  interceptor  missile  flight  path,  line- 
of-sight  angles  and  their  derivatives.  Estimation  of  the 
target  position,  velocity  and  acceleration  is  'ot 
required.  We  have  validated  the  effectiveness  of  the 
designed  sliding  mode  controllers  through  numerical 
examples. 

Introduction 

Kinetic  energy  kill  vehicles  are  employed  as 
interceptors  in  theater  missile  defense  (TMD)1.  A 
KEKV  destroys  an  incoming  ballistic  missile  as  they 
impact.  The  interceptor’s  guidance  and  control  system 
must  provide  autonomous  target  tracking  and  homing 
maneuvers  to  achieve  a  very  accurate  intercept.  A  two- 
dimensional  missile-target  engagement  geometry2  is 
shown  in  fig.  1,  where  XOZ  is  the  inertial  frame,  V M 
and  Vr  are  velocity  magnitudes  of  the  missile  and  the 
target  (m/s),  CGM  and  CGr  are  centers  of  gravity  of 
the  missile  and  the  target,  X  is  the  line-of-sight  angle 
(rad),  y  and  p  are  flight  path  angles  of  the  missile 

and  the  target  respectively  (rad),  Rm  is  a  length  of  the 


line-of-sight  or  range  (m),  a  is  angle  of  attack  (rad), 
nc  is  the  acceleration  missile  command  (mis2),  nT  is 

the  target  acceleration  ( m/s2 ),  ud  and  ua  are  divert 
and  attitude  control  forces  respectively. 


Fig.  1  Two-dimensional  missile-target  engagement 
geometry 

The  jamming  signal  is  a  random  noise  generated  by  the 
target.  Traditionally,  the  acceleration  missile  command 
nc  is  calculated  as  proportional  or  augmented 
proportional2  navigation  laws,  advanced  guidance 
law2,3,  linear  quadratic  optimal  guidance  law4,  etc.  The 
acceleration  missile  command  nc  is  then  executed  by  a 
flight  control  system  to  minimize  the  miss  distance  (the 
distance  between  the  target  and  interceptor). 

The  simplest  proportional  navigation  law2  can 
be  expressed  in  a  format 

nc  =  NVclX  ,  (1) 

where  N  e  [3,5]  is  the  navigation  ratio,  and  VCL  is  a 
magnitude  of  the  missile-target  closing  velocity  (m/s). 
It  is  well  known2  that  implementing  the  proportional 


Distribution  Statement  A.  Approved  for 
public  release;  distribution  is  unlimited. 


1 

Unclassified 


DTIC  QUALITY  INSPECTED  4 


Unclassified 


navigation  law,  Eq.  (1),  requires  the  seeker  to  estimate 
X  as  well  as  VT ,  p ,  nT  .  To  estimate  the  values  of  the 
target  parameters,  high  order  Kalman  filters2,3  are 
usually  designed.  These  filters  should  also  address 
nonlinearities  of  the  missile-target  kinematics  and 
dynamics.  A  maneuvering  target  (for  instance, 
submunition)  complicates  predicting  target  motion. 
Sometimes  these  estimations  are  done  by  postulating 
the  sinusoidal  target  weave  motion5.  A  hostile 
environment  affects  the  interceptor’s  motion.  In 
particular,  unknown  disturbances  (for  instance,  wind 
torque  during  atmospheric  flight)  and  parametric 
uncertainties  including  control  system  failures  (for 
instance,  damage  of  the  aerodynamic  surfaces)  affect 
the  interceptor's  motion  and  tracking  dynamics.  A 
jamming  signal  (noise)  which  is  generated  by  the  target 
can  create  difficulties  for  the  seeker  in  estimating  the 

values  of  X  ,  VT ,  p  and  As  a  result,  the 
guidance  and  control  system  of  the  KEKV  becomes 
cumbersome  and  extremely  sensitive  to  the  accuracy  of 
parameter  estimation  for  target  motion  and  unknown 
disturbances. 

One  of  the  promising  nonlinear  control 
algorithms  that  is  robust  to  disturbances  and  plant 
uncertainties  is  Sliding  Mode  Control6  8.  The  purpose  of 
a  Sliding  Mode  Controller  (SMC)  is  to  drive  a 
nonlinear  plant's  trajectory  onto  a  prescribed  (user- 
chosen)  surface,  named  the  sliding  surface,  in  the  state 
space  and  to  maintain  the  plant’s  state  trajectory  on  this 
surface  thereafter.  The  motion  of  the  system  on  the 
sliding  surface  is  called  the  sliding  mode.  The  equation 
for  the  sliding  surface  must  be  selected  such  that  the 
system  will  exhibit  the  desired  (given)  behavior  in  the 
sliding  mode  independent  of  unwanted  parameters 
(plant  uncertainties  and  disturbances).  Usually  the  SMC 
is  characterized  by  a  high  frequency  switching  control 
function  that  provides  the  system's  motion  in  the  sliding 
surface6"8.  The  discontinuous  nature  of  a  SMC  is 
suitable  KEKV  control  systems  that  use  divert  and 
attitude  thrusters.  A  SMC  can  be  implemented  in  a 
pulse-width  modulation  (PWM)  format  to  enforce  a 
certain  frequency  of  thruster  firing.  However,  guiding 
the  KEKV  by  aerodynamic  surfaces  requires 
continuous  implementation  of  a  discontinuous  SMC9,10. 
SMCs  have  been  successfully  applied  to  address 
significant  plant  uncertainties  and  disturbances  for 
various  space  and  flight  control  problems’11,12. 

SMC  Design  for  KEKV 

Robust,  simple  guidance  and  control 
algorithms  for  kinetic  energy  kill  vehicles  operating  in 
hostile  environments  can  be  developed  in  sliding 
modes.  The  simplicity  of  the  SMC’s  design  is  based  on 


the  fact  that  estimating  the  target’s  parameters  is  not 
required.  Only  X ,  y ,  X  and  y  must  be  measured  or 
estimated  for  the  interceptor.  The  robustness  of  the 
SMCs  has  been  demonstrated  on  system  operation  in 
sliding  mode6"11.  In  designing  the  SMC  for  a  KEKV, 
two  goals  are  proposed: 

1.  Direct  the  velocity  vector  of  interceptor,  V*.  to 
the  target  by  following  the  line-of-sight  angle  X  to 
the  interceptor  flight  path  angle  y  regardless  of 
disturbances  and  uncertainties, 

2.  Reduce  drag  force  by  directing  the  interceptor  axis 
of  symmetry  to  the  center  of  gravity  of  the  target 
regardless  of  disturbances  and  uncertainties. 

We  have  evaluated  SMC  designs  for  two  types  of 
KEKVs,  a  ballistic  missile  and  a  cruise  missile  type. 

SMC  Design  for  a  Ballistic  Type  KEKV 

A  simplified  mathematical  model  is  obtained 
for  the  ballistic  type  interceptor  in  the  pitch  plane. 

Mathematical  Model  of  a  Plane  Motion 
A  simplified  translational  motion  of  the 
ballistic  type  interceptor  in  the  pitch  plane  (fig.  1)  is 
described  as  follows: 

=-—(ua  +Wrf)sin0  +  — , 

m  m  (2) 

vz  =  —  (ua  +ud)cos0-g  +  ^-, 
h  m  m 

where 

Vv=K.v,r.  m  is  the  mass  of  the  interceptor  (kg), 
'V  =  {ysx,  y/z  }r  is  the  vector  of  disturbance  and 

aerodynamic  forces  ( N ),  and  g  -  9.81m  Is2 . 

A  simplified  rotational  motion  of  the  ballistic 
type  interceptor  in  the  pitch  plane  (fig.  1)  is  described 
as  follows: 

I #  =  q> 


where 

q  the  a  pitch  rate  (rad/s),  J  is  the  moment  of  inertia  of 
the  interceptor  ( kg-m 2),  r  is  the  distance  between 
application  of  the  attitude  control  force  ua  and  the 
interceptor  center  of  gravity  CGM  ,  and  y/^  is  the  sum 
of  aerodynamic  and  disturbance  torque. 

The  output  equations  are 
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yx  =  y  =  tan 

yi  =  °- 


(4) 


Problem  Formulation 

To  address  the  first  goal  for  the  SMC,  we  will 
provide  the  following  tracking  motion: 

\\m\yu(t)-yx(t)\  =  \im\X(t)-y(t)\  =  0,  (5) 

f— >00  t— >°0 

where  d  denotes  the  desired  output. 

To  address  the  second  goal  (attitude  control 
problem),  the  SMC  is  designed  to  fit  a  desired  tracking 
motion 

lim| y2d(t)  -  y2(t)\  =  lirn|A^  -  6(t)\  =  0  . '  (6) 

We  can  observe  that  the  rotational  motion  in  not 
effected  by  divert  control  force  ud .  Therefore,  the 
second  goal,  Eq.  (6),  will  be  addressed  first  via  the 
design  of  an  attitude  SMC,  ua  .  The  first  goal,  Eq.  (5), 
will  be  addressed  next  through  the  design  of  a  divert 
SMC,  ud . 


SMC  Design 

To  address  the  attitude  control  problem,  Eq. 
(6),  the  following  sliding  surface  was  designed 

a  2=e2+c2e2=0,  (7) 

where  e2  =X(t)~0(t) . 

The  value  of  a  parameter  c2  will  be  identified 
to  provide  the  desired  eigenvalue  to  a  homogeneous 
linear  time-invariant  differential  equation  (7)  that 
describes  the  tracking  dynamics  of  a  rotational  motion 
on  the  sliding  surface. 

The  control  law  ua  is  designed  to  provide  a 
finite-time  convergence  of  the  system’s  Eq.  (3) 
trajectory  to  the  origin  in  the  cr2  -  subspace.  The 
system's  Eq.  (3)  motion  in  the  cr2  -  subspace  is 
derived  as  follows: 


(8) 


&2  —  X(t)  +  c2e2  j  “*■  j  ua' 

A  candidate  to  the  Laypuniv  function  is 
formed  as  follows: 


CTo 

and  its  derivative  is  identified 
02  =  *2*2  =  ^2 


(9) 


^ ..  w  r  ^ 

Mf)  +  c2e2-*±+jua 


(10) 


v  ^  J 

Providing  a  finite-time  convergence  to  a2  =  0  the 

control  law  ua  is  designed  to  meet  the  following 
inequality6'8: 


(11) 

is  identified  as 


Q2=a2a2<-p\a2\,  p>  0 
Equivalent  attitude  control  ua^ 
follows6'8: 


Uaeq  +"y 

Meeting  inequality  (11)  and  assuming 
.  ~ua  <N>0, 

eq  ueq 

the  following  SMC  is  designed: 

K  =“am  ~P2signa2, 


(12) 

(13) 

(14) 


where  u .  is  u„  estimate,  and  p2  >N  +  —  p. 

ueq  ueq  j. 

Directing  the  interceptor  velocity  vector,  VM, 
to  the  target,  the  divert  SMC,  ud  ,  is  designed  to  meet 
Eq.  (5).  Differentiating  yx  we  obtain 


cos  a  /1C. 

yt  =Q(.)+  .  ,ud,  (15) 

I 

where 


vw|  =  Vv^  +vz  and  a  =  6-y 


The  following  assumptions  are  made 

1 .  ina  reasonable  flight  domain, 

m\yM\ 

2.  corresponding  internal  (zero)  dynamics  are 
stable13. 

Then  the  relative  degree  is  equal  to  one,  and  the  sliding 
surface  is  designed  as  follows: 

o-j  =  gj  =  0 ,  (16) 

where  ex  =  X(t)-y(t)  . 

It  is  obvious  that  after  the  sliding  surface  in 
Eq.  (16)  has  been  reached,  the  conditions  of  Eq.  (5)  will 
be  met.  The  control  law  ud  is  designed  to  provide  a 
finite-time  convergence  of  the  system's  Eq.  (2) 
trajectory  to  the  origin  in  the  crx-  subspace.  The 
system's  Eq.  (2)  motion  in  the  ax  -  subspace  is  derived 
as  follows: 


cos  a 


(17) 


A  candidate  to  the  Laypuniv  function  is 
formed  as  follows: 

cri2 


Qi  =_^_>0, 


(18) 
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and  its  derivative  is  identified 


Qi  =  = 


Providing  a  finite-time  convergence  to  <jx  =  0 ,  the 
control  law  ud  must  be  designed  to  meet  following 
inequality6'8: 

Qi  =  o-jd-,  <  -p\<J\  |,  p>  o  (20) 

The  equivalent  divert  control  ud^  is  identified6'8as 


=  ~Ua  + 


m 

cos  a 


y/x  . 

cos^  +  -^-sin^  + 
m 


(21) 


Assuming 

ud  -ud  |</5>0,  (22) 

ueq  ueq  J 

the  following  SMC  is  designed  to  meet  the  inequality 
described  in  Eq.  (20): 

ud  =  ^  +  PiSign&i  >  (23) 


where  ua  is  ua  estimate,  and  px  >P  + 


m\y 


M\ 


cos  a 


Simulation 

The  following  numerical  values  represent  a 
generic  ballistic  missile  KEKV: 

m-lOkg ,  J  =  I0kg-m2 ,  r  =  1  m.  (24) 

The  maximum  values  of  attitude  and  divert  thrust  are: 


Uamm  =  500 N,  Udmm  =  40000  N .  (25) 

The  following  interceptor  initial  conditions 
were  used  during  simulations: 

|z(0)=  50,000m,  x(0)  =  10,000m,  0(0)  =  -0.5 rad,  ^ 
s  (26) 

[®(0)  =  0,  F,(0)  =  2,000 m/s,  Fz(0)  = -1,000 m/s 

The  initial  position  and  velocity  of  a 
maneuvering  target  (submunition)  were  chosen  as: 


J  zT  (0)  =  40,000  m,  xT  (0)  =  35,000  m, 

\vxT (0)  =  -750 m / 5,  VzT (0)  =  -100 mis  " ' 1 

The  attitude  SMC,  ua ,  was  designed  in 
accordance  with  Eqs.  (7)  and  (14).  This  is 
ua  =  500 sign<r2  5  cr2  =  e2  +  2e2  =  0  .  (28) 

The  divert  SMC,  wj,  was  designed  in 
accordance  with  Eqs.  (16)  and  (23).  This  is 
ud  =40,000<?/g?2o-1,  <jx  -ex  (29) 

The  results  of  the  system  "interceptor-targef 
simulation  without  disturbances  and  sensor  noise  are 


presented  in  figures  2  through  8.  Very  accurate  kill  and 
attitude  alignment  is  demonstrated.  The  results  of  the 
systeiris  simulation  with  disturbances,  including 
jamming  signals,  and  censoring  noise  are  presented  in 
figures  9  through  11.  The  line-of-sight  angle,  A,  is 
accurately  followed  by  the  interceptor  flight  path  angle, 
y ,  and  the  pitch  angle,  6 ,  regardless  of  the 
disturbances  or  measurement  noise. 


Angles  (rad) 


Fig.  2  Tracking  of  line-of-sight  angle  without 
disturbances  and  noise  of  measurement 

Length  offtie-or-sight  (meters) 


Fig.  3  Length  of  line  of  sight  without  disturbances 
and  noise  of  measurement 
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SMC  Design  for  a  Cruise  Type  KEKV 

We  now  consider  a  strategic  air  defense 
missile  interceptor  (cruise  type),  using  tailfln  control. 


Mathematical  Model  of  a  Plane  Motion 

The  following  simplified  mathematical  model 
represents  the  pitch-plane  angular  equations  of  motion 
of  the  cruise  type  interceptor14. 


6 

k 

V 


0 


1 


qc-S 
J  ^ 


■K|  * 

0 


0 

qc-S 

~T  ^ 


0 

0 

0 

+ 

q-l-S 

J  Cn* 

cog' 

n.K 

1 

5 

_ 1 

\e' 

*=[ o  o  1 \q 


(30) 


where  S  is  the  tailfin  cross  section  area  (m2),  q  is  the 


dynamic  pressure  (N /  m2 ),  c  is  the  missile  diameter 
(m),  Cma,  Cmq,  Cmq,  CZa,  C2s  are  aerodynamic 
coefficients,  m  is  the  mass  of  the  KEKV  (kg). 


Problem  Formulation 

Having  only  one  control  input:  deflection  8 
of  a  tailfin,  we  can  address  only  one  goal:  following  the 
line-of-sight  angle,  X ,  by  the  interceptor  flight  path 
angle  y.  Addressing  this  goal  we  will  provide  the 
following  tracking  motion: 

\\m\yld(t)-y^(t)\  =  \im\X(t)-y(t)\  =  0 

t — >00  t~~ >°0 

by  a  corresponding  design  of  a  control  input:  deflection 
5  of  a  tailfin.  This  design  will  be  accomplished  in  a 
continuous  SMC  format9,15. 

We  can  observe  that  the  control  input  8 
enters  two  equations  of  the  system  in  Eq.  (30).  This  fact 
can  lead  to  a  nonminimum  phase  case13.  It  is  well 
known  that  traditional  sliding  mode  control  is  not 
applicable  to  nonminimum  phase  systems6,7.  Sliding 
mode  control  that  is  applicable  to  nonminimum  phase 
systems  is  developed  on  the  basis  of  dynamic  sliding 
surface  design16. 

Basics  of  Nonminimum  Phase  Output 

Tracking  in  Dynamic  Sliding  Surfaces16 


The  following  nonlinear  nonminimum  phase 
plant  is  considered 

x  =Ax+F(x,t)+bu,  y=Gx,  (31) 

where  x  e  iHn  is  the  state  vector,  u  e  SH1  is  the  control 

function,  yeM1  is  the  controlled  output,  A,b,G  are 

the  constant  matrices  of  corresponding  dimensions, 

{. A,b }  is  a  controllable  pair,  F(x,t)e  is  a  nonlinear 

time-dependent  vector-function,  so  that 

F(x,t)  =  Fx(x)+F2(t) ,  where  Fx(x) :  ||JF1fx)|[  <  Nx  ||*|[  is  a 

matched  nonlinear  function  and  F2(t) : 

\F2/t)\<N2JVi  =  ln  is  a  smooth  enough  vector- 

disturbance  with  the  positive  constants  NlfN2,i  * 

We  wish  to  specify  the  sliding  mode  controller 

\u+ ,3(x,e,t)>  0, 
u=\ 

[tt  ,3(x,eJ)<0, 
where 

3(x,e,t)=  0  (33) 

is  an  equation  of  sliding  surface  as  a  dynamic  operator 
acting  on  state  variable  vector  x  and  output  tracking 

error  e(t)=y* (t) - y(t)  ;  u+ ,u~  -  continuous  functions  of 
x,t ,  in  order  to  accomplish  the  following  goal: 

•  the  output  y(t)  of  the  nonlinear 
nonminimum  phase  system,  Eq.  (31),  must 
asymptotically  track  the  given  reference  output 

trajectory  y*(t)  in  a  sliding  mode.  That  is 
lim|ly*^-y(7J=lim||e(7j:z::0 ,  and  the  output  tracking 

f->oo!l  II  f— >00"  11 

must  be  linear  with  given  eigenvalues  placement  in  the 
dynamic  sliding  surface  from  Eq.  (33). 

Utilizing  a  nonsingular  transformation6,7 


(32) 


A/1 


M 


\k,  M- 


,  MB= 


AnA2 


A2lA22 


b2*0,  MF= 


MAM1  = 

gm -1  =[g’:g2] 


Ml 
M2 

m 

A A(x\x2)+f2(t) 
of  the  system  in  Eq.  (31)  to  a  regular  form6,7,  the 
system,  Eq.  (3 1),  is  rewritten  in  the  following  format: 

xl  =  Auxl  +  Anx2  +  fx(t) 

x2  =A21  x1  +A22x2  +  lsA(xx  ,x2)  +  f2(t)  +  b2u  (34) 

y  =  Gxxx  +G2x2, 

where  x1  eRn~l,x2  eRl,b2  ^0.  Considering  a  state 
variable  x2  as  a  virtual  control,  the  transfer  function 
Y(s) 


X2(s) 


=  Gj  [si  —  Ax  |  ]  1  Al2  +  G2 


(35) 
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is  assumed  to  have  zeros  in  the  right  hand  half  of  the 
complex  plane.  It  is  obvious  that  a  direct  application  of 
a  conventional  sliding  mode  controller  to  an  output 
tracking  in  the  system  (34)  is  not  possible  because  of 
the  nonminimum  phase  nature  of  the  plant6,7,13 

Considering  a  state  variable  x 2  as  a  virtual 
control  input  at  the  first  step  of  the  sliding  mode 
controller  design,  the  dynamic  sliding  manifold  is 
introduced  as  follows: 

3  =  x2  +cr  =  0  .  (36) 

The  function  cr  is  designed  as  a  dynamic  operator, 
acting  on  the  output  tracking  error  e  : 

a  =  W(s)e  =  %^e,  (37) 

P(s) 

where  Q(s)  and  P(s)  are  polynomials  of  the  Laplace 
variable  s . 

Assuming  that  the  sliding  mode  exists  in  the 
system,  Eq.  (34),  on  the  dynamic  sliding  surface,  Eq. 
(36),  the  motion  of  the  system  in  Eq.  (34)  in  this 
surface  is  described  as  follows6,7: 

i1  =Anxl  -  Ancr  +  fx(t),  y  =  Glxl-G2cr  (38) 
The  evolution  of  e(t)  in  the  system,  Eq.  (38), 
was  identified  in  operator  notation.  This  is 
h)  -  G2Q(s)  -  G,  (5/  -  An)-1  AnQ{s)\{s)  = 

=  P(s)Y*  (s)  -  G,  (si  -  An )~:  P(s)Ft  (s). 

The  polynomials  P(s)  and  Q(s )  are  designed 
in  order  to 

•  provide  any  desired  roots’  placement  to  the 
characteristic  equation 

P(s)  -  G2Q(s)  -  G,  {si  -  A, , )"'  A12Q(s)  =  0 ,  (40) 

•  reject  the  effect  of  the  reference  profile  y*(t)  and 
the  unmatched  disturbance  fx(t)  to  the  output 
tracking  error  e  in  a  steady  state. 

The  rejection  condition  is  derived  as  follows: 

=ijng  kzAi  rgfww  =o.  (4i) 

*-°  P(s)-G2as)-Gx{sI-4i)  42&S) 

A  control  function  u  is  designed  to  meet  the 
well  known6,7  sliding  mode  existing  condition 

3-3<-p|3|,  p  >  0  (42) 

Eq.  (42)  is  calculated  for  the  system  in  Eq.  (30).  This  is 

u+  < — -[<j+ p+An^  +  A22x2+AA(x\x2)  + /2(t)] 

bl  (43) 

u  >  — — Lt- p+A2Ix'  +  A22x2  +AA(x\x2)+ /2(/)] 
bi 

where  b2  is  assumed  to  be  positive.  To  make  Eq.  (43) 
realizable,  we  need  to  make  cr  bounded.  Obviously 


j  As  Pm(s)\  1 

1  -  G2m  -  g{  {si-ax  Al2a4 

is  bounded  if 

•  the  output  tracking  profile  y*(t)  and  the 
disturbance  fx(t)  are  bounded, 

•  the  denominator  of  Eq.  (44)  does  not  have  roots  in 
the  right  half  of  the  complex  plane,  and  has  at  most 
one  root  in  the  origin. 

Consequently,  the  sliding  mode  exists  in  the 
nonminimum  phase  system  (31)  via  SMC  Eqs.  (32)  and 
(43)  in  the  dynamic  sliding  manifold 

3  =  M2x  +  Q&-e  =  0,  (45) 

P(s) 


SMC  Design 

The  mathematical  model,  Eq.  (30),  was 
considered  for  a  generic  strategic  air  defense  missile 
interceptor  at  Mach  4.5  and  altitude  10  km.  This  is 

e=q 

<7=-3873(9+3874y-2.5^-4382<5+ f2(f) 
y=l-2390-1.239!r-O.OO68cosy+O.471(5+^(O’  1  ' 


l y=r 

This  system  is  transformed  to  a  regular  form 
of  Eq.  (34)  as  follows: 

2  =  -0.823z + 0.823(7  -  0.001 9q  - 
-0.0068cos£  -0.001075/)  +/,  (0  -0.00107^(0 
■0  =  q  (47) 

q  =  3873z  =  3S14d-0.675q-43&2S+f2(t) 


y  =  z-0.00l075q 

The  transfer  function  Eq.  (35)  identified  for 
the  system  in  Eq.  (47)  is 


^  = -0.001075^ 


-26.5)(s  +  28.9) 


(48) 

q(s)  5(5  +  0.823 

which  has  a  zero  in  the  right  hand  side  of  the  complex 
plane  and,  therefore,  is  nonminimum  phase.  The 
dynamic  sliding  surface,  Eq.  (36),  is  designed  as 
follows: 


3  =  #  +  <x  =  0.  (49) 

The  system's  Eq.  (47)  motion  in  the  sliding  surface  of 
Eq.  (49)  is  identified  in  accordance  with  Eq.  (38).  This 
is 


z  =  -0.823 z + 0.823 6 + 0.00 1 9cr + ^z,  cr,  t) 

<e  =  -a  (50) 

y  =  z -0.00 1075# 


where  ^(z,^,/)  =  -0.0068cos^+ fx (t)~ 0.001 07 5f2 (t)  . 
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Assuming  y*(t)  =  yld  and  y/{z,q,t)  are 
piece-wise  constants,  the  following  polynomials  P(s ) 
and  Q(s)  were  identified 


(51) 


f  Q(s)  =  -287. 19s 2  -  2259.365'  -  8399.39 
[P(s)  =  s(s  + 16.91) 

It  is  easy  to  show  that  the  identified  P(s)  and 
Q(s )  meet  Eq.  (41),  and  make  the  characteristic 
equation  of  the  system  Eq.  (50) 

[POr)  +  O.OO19750(s)>2+  (52) 

+  [0.823/,(5)  +  O.OO2750(j)]- 0.8230(5)  =  0 
to  satisfy  ITAE  criterion 

s4  +2Acos3  +3.4  co2s2  +2.7  a)3s  +  co4  =0  (53) 

with  co  =  10  rad  I  s  .  So,  the  following  dynamic  sliding 
surface  is  designed 

^  287.19s2  +  2259.36^  +  8399.39  _ 

3  =  q - e.  (54) 

.  s(s  +  16.91) 

The  control  law  for  Eqs.  (32)  and  (43)  is 
designed  in  a  simplified  format.  This  is 


8  =  SmzxsignZ , 
where 


^max  ^  1 8 eq  |  P 


'eq 


P>  0, 

=  a  +  387.3z  -387.30  -  0.675?. 

Eliminating  control  chattering, 


(55) 


(56) 


the 


discontinuous  control  law  in  Eq.  (55)  is  substituted  by 
3 


8  =  8msat- 


£  >0. 


(57) 


Simulation 

The  cruise  missile  type  interceptor  Eq.  (46) 
was  simulated  with  the  dynamic  SMC  Eqs.  (54)  and 
(57)  and  £max  =  0.532  rad  and  s  =  1 .  A  low  pass  filter 

with  a  transfer  function  - - -  filters  the  dynamic 

1  + 0.03s 

sliding  surface.  The  results  are  shown  in  figures  12 
through  14.  It  is  observed  that  the  tracking  of  the  sensor 
with  noise  line-of-sight  angle  profile  by  the  flight  path 
angle  is  insensitive  to  the  disturbances  f(t)  and 

flit). 

Conclusion 

We  have  shown  the  development  of  a  high 
performance  hit-to-kill  interceptor  control  system, 
addressing  hostile  environment,  which  includes  noises 
and  disturbances,  based  on  sliding  mode  control 
techniques.  We  have  addressed  the  control  of  both 
ballistic  and  cruise  missile  type  interceptors.  After  the 
interceptor  flight  path  angle  has  captured  the  line-of- 


sight  angle,  the  SMC  will  guide  the  interceptor  along 
the  line-of-site  angle  until  the  target  missile  is 
impacted.  The  SMC  control  algorithms  we  have 
developed  do  not  require  estimations  of  the  target 
position,  velocity  and  acceleration.  The  result  of  this 
development  is  robust,  relatively  simple  controllers  for 
TMD  missile  applications.  Future  work  needs  to  be 
accomplished  to  develop  high  fidelity  SMC  algorithms 
for  incorporation  into  TMD  missile  simulations  and 
systems. 
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Output  Tracking  (rad) 


Fig.  12  Tracking  of  line-of-sight  angle  with 
disturbances  and  noise  of  measurement 
(cruise  missile  type  interceptor) 


Fig.  13  Control  input:  tailfin  deflection  with 
disturbances  and  noise  of  measurement 
(cruise  missile  type  interceptor) 


Disturbances  (nd^sec) 


12'  1 — 1 _ i _ i _ i _ i _ i _ i _ i _ i _ i - ! 

0X23456789  10  1112 

Time  (sec) 


Fig.  14  Disturbances  (cruise  missile  type  interceptor) 
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